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CHfford quantum cellular automata (CQCAs) are a special kind of quantum cellular automata 
(QCAs) that incorporate Clifford group operations for the time evolution. Despite being classically 
simulable, they can be used as basic building blocks for universal quantum computation. This is 
due to the connection to translation-invariant stabilizer states and their entanglement properties. 
We will give a self-contained introduction to CQCAs and investigate the generation of entanglement 
under CQCA action. Futhermore, we will discuss finite configurations and applications of CQCAs. 



I. INTRODUCTION 

In this work we will give an introduction to ClifFord 
quantum cellular automata (CQCAs) [J [3| and study 
their entanglement generation properties. Quantum cel- 
lular automata (QCAs) Q are a quantum computational 
model that only requires global control for the time evo- 
lution. The only part of a QCA computation that needs 
local control is the preparation of the input state, which 
in the case of a QCA capable of universal quantum com- 
putation 0,01 encodes the program and the data input. 
But not all QCAs are universal quantum computers. The 
CQCAs we deal with here are even classically simulable. 
Yet they are of great use in several quantum computa- 
tional schemes. Furthermore, the classical structure al- 
lows an in-depth analysis of their time evolution, which 
can give us a good intuition about the properties of gen- 
eral QCAs. Most of the results presented here are also 
included in In this work we focus on entanglement 
generation and try to present the theory of CQCAs in a 
self-contained manner. For a more detailed and mathe- 
matically concise formulation we refer to [l| and . 

We will begin with a short introduction to QCAs and 
CQCAs as well as the classical description of CQCAs. 
Then we will introduce translation-invariant pure sta- 
bilizer states and study their entanglement in a bipar- 
tite setting. Putting the results of the first two sections 
together, we will investigate the evolution of stabilizer 
states under CQCA action and derive a simple expression 
for the generation of entanglement by CQCAs acting on 
stabilizer states. Namely, the entanglement generation is 
linear and depends only on the degree of the trace of the 
polynomial matrix representing the CQCA. We will then 
make some short remarks on CQCAs over finite chains 
and finally give a short introducti_on to applications of 
CQCAs in quantum computing 



II. 



CLIFFORD QUANTUM CELLULAR 
AUTOMATA 



on a translation-invariant lattice of quantum systems. 
So, in each discrete step of the time evolution of our 
lattice of quantum systems, we will do the same thing 
and we will do it in a way independent of the position 
in our lattice. Because we deal with infinite lattices, e.g. 
spin chains, it is convenient to look at the time evolution 
of observables rather than states to circumvent mathe- 
matical problems with Hilbert spaces on infinitely many 
tensor factors. Using the observable evolution, we get 
the nice property that a QCA is uniquely determined by 
the image of the observables on a single site Q. We call 
this local operation T; if it acts on the observables of the 
ith system. As we defined our QCA to be translation 
invariant, all T^ are the same and we pick one, namely 
To, as a representative. We now introduce another prop- 
erty of QCAs: locality. In each time step we only allow a 
finite propagation of information on the lattice, i.e. two 
states that only differ on one site may only differ on a 
neighborhood of finitely many sites around the original 
site after the application of the QCA. This means that 
the image of the observables at site i is a subset of the the 
observables of a finite number of sites surrounding site i. 
We call this set the neighborhood N of the automaton. 
This is depicted in Figure [1] 



A. Reversible Quantum Cellular Automata 

A reversible quantum cellular automaton (QCA) T is 
a translation-invariant reversible discrete time operation 



In mathematical terms we say that a QCA is a lo- 
cal, translation-invariant automorphism of the algebra 
of observables 21. An automorphism T has the prop- 
erty T{AB) = T(^)T(B), VA,S e 21. The local op- 
eration To is a homomorphism from 2li to 2l(A/' -I- «), 
where J\f + i denotes the neighborhood of the ith cell 
and 2ti the observable algebra of the cell. A homomor- 
phism preserves commutation relations in the sense that 
[T{A),T{B)] = T[A,B] holds. But in our local rule To, 
information about the evolution of neighboring cells does 
not play a role. However, for T to be an automorphism, 
the images of two commuting observables still have to 
commute after the evolution. So, we impose this condi- 
tion on To by demanding that all observables in T(2ti) 
commute with all observables in T(2ti-f2-), x 0, where 
T(2li) denotes the algebra of the images of observables 
localized on cell i. 



To : M ^ M + A/" 



□□□□□□□□□□□□□□□□ ' + ' 



Figure 1: Time evolution of a CQCA. The images of the observables on system s^i are contained in the set of observables of 
the systems ^ + A/" in the next time step. 



B. Clifford Quantum Cellular Automata 

Now we restrict ourselves to a speeial kind of QCAs, 
the Clifford quantum cellular automata (CQCAs), which 
employ the Clifford group operations. Clifford operations 
map tensor products of Pauli matrices to tensor products 
of Pauli matrices times a phase. Thus, the local rule Tq 
of a CQCA maps single Pauli matrices to tensor products 
of Pauli matrices times a phase. Furthermore, we only 
consider one-dimensional lattices, i.e. spin chains. We 
will illustrate this with an example. 

Example II. 1. The so-called glider CQCA Tq has the 
following local rule: 



Zi. 

•X, 



where Xi and Zi denote the Pauli matrices at system i. 
The image of Yi follows from the product of the images of 
Xi and Zi, because we require Tq to be a homomorphism: 

TGfim^-Z^-l<^^Y,®Z,+ l. 

The global transformation Tq has to be an automor- 
phism. We already constructed the image of Yi accord- 
ingly. Now we have to check if the commutation and 
anti- commutation relations between neighboring cells are 
preserved. We have 

[Tg,oX\Tg.oX'] = [Z\Z^]^0, 



[Tg.oZ\Tg,oZ^] 

[Z'~^ (g) X' (g) Z''+\Z^~^ > 

0, 



. X^ (g) Z^+'] 



[Tg,oZ\Tg^oX^] = [Z'"'( 
{TgmZ\Tg,oX'} = {Z'-' 



<X' 
iX' 



Z''+\Z^]^Q,z^3, 
iZ'+\Z'} = 0, 



so our local rule extends to a global automorphism. As 
the global rule acts in exactly the same way on single-site 
Pauli matrices as the global rule, we will from now on 
always use the global rule Tq. 

By neglecting a global phase, we can also think of the 
CQCA Tg as a classical automaton acting on the la- 
bels (1 = X, 2 = y, 3 = Z, = 1) of the Pauli matrices. We 



define the operation Q, which has the following prop- 
erties to resemble the multiplication of Pauli matrices: 
101 = 0, iQ j = k for i,j,k = 1,2,3 with i ^ j ^ k and 
O0i = i0O = i for i = 0, . . . , 3. Now we can illustrate 
the evolution of the observable Z~^ ®Y^ ® X^ as follows 
(the underlined labels are situated at the origin.): 



Tg(321) = Tg(3OO)0Tg(O2O)0Tg(OO1) 

3 1 3 

_ 3 2 3 

~ ^ 3_ 

= 3210 



(3 21) 



We can see that this observable only moves on the lattice 
by one step. We call observables with this property "glid- 
ers" and the automata which allow for such observables 
glider automata. Figure [H gives an explanation for this 
behavior. 



1. Classical Description of CQCA 

It is well known that Clifford operations are classically 
simulable in an efhcient way Q . We use this to derive a 
classical description of CQCAs, which greatly simplifies 
our further analysis. Heuristically speaking we replace 
the CQCAs acting on observables by a classical cellular 
automaton (CA) acting on the labels of Pauli matrices 
as we did with the example CQCA in the last section. 
Mathematically we do the following: first we introduce a 
function w that maps vectors over Z2 to Pauli matrices^: 

a: = w(i,o), r = iw(i,i), z = w(o,i), i = w{q,q). 

(1) 

This way we gave the Pauli matrices classical labels from 
a vector space. A quantum cellular automaton mapping 
Pauli matrices to Pauli matrices can now be described 
as a classical automaton acting on the labels. But, be- 
fore we introduce the classical description of CQCAs, we 
have to introduce tensor products of Pauli matrices to 
our classical description. To do that we extend the func- 
tion w in a way that it maps vectors of binary strings to 



is the finite field consisting of and 1 with addition modulo 



Figure 2: Time evolution of the glider CQCA Tq'- On the left side the individual evolutions of single Z and X observables are 
shown. They both generate "checkerboards" which are shifted against each other by one time-step. By shifting the original X 
also in space we overlap the checkerboards such that they cancel out leaving only a small strip in the space-time picture. This 
strip is generated by the original observable moving on the spin chain. Such observables are called gliders. 



tensor products of Pauli matrices: 



w(0 



The ones in account for the PauU X matrices, the ones 
in for the Y matrices, and if £,xixo) = £,y{xo) = 1, 
there is a Y at xq: 



110 10 
110 



■ t<^X<^iY<SiZ<SiX<^t 



The binary strings arc infinitely long, but, as wc always 
deal with localized observables, only finitely many posi- 
tions have non-zero entries. We omit the zeros before and 
after the observable to get a more convenient notation. 
Similarly, we will omit the identities before and after the 
observable. The binary strings form a vector space over 
Z2, which we call phase space. 

The addition in the classical description is commuta- 
tive, so we have to encode the commutation relations in 
an auxiliary function. Our operators fulfill the condition 

w(e + r/) = (-l)''-«^w(0w(7^), 

which is known from the Pauli matrices as 

(^k = ieijkcncTj, 

where ai = X, <J2 = Y , (T3 = Z, and e^fc = 1 if k) is 
an even permutation of (1, 2, 3), Sijk = —1 if («, J, k) is an 
odd permutation of (1,2,3), and Sijk = otherwise. In 
our definition we have no i, because we chose w(l, 1) = 
lY . Therefore, the w(^) fulfill the commutation relations 

= (-l)«^''^~«^''-'''w(?7)w(0. (3) 

In both cases terms of the type S.xVz SuTC scalar products 
where the addition is carried out modulo 2.^ We call the 
function 



<^{tv) = ^xVz - ^zVx 



(4) 



the symplectic form. Thus, the commutation relations 
are encoded in the symplectic form. 

To further simplify our classical description, we trans- 
form the binary strings to Laurent polynomials using 



(5) 



We use these polynomials as abstract objects without 
ever evaluating them for some value of u. Our example 
observable then looks like 



110 10 
110 



+ 1 + 
l + u 



We change the function w in a way that it maps vectors 
of Laurent polynomials to tensor products of Pauli ma- 
trices. By a slight abuse of notation, we will also name 
it w. 

Now wc want to describe CQCAs in this picture. A 
CQCA T maps every Pauli matrix to a tensor product of 
those. Because it is an automorphism and determined by 
the local transformation Tg, we can calculate the image 
of any observable from the image of the Pauli matrices it 
can be decomposed into. Of course, the images of neigh- 
boring Pauli matrices in general overlap and have to be 
multiplied. Because we required them to commute, this 
multiplication is uniquely defined. In the phase-space 
picture with binary strings this multiplication becomes 
an addition. The CQCA's action is thus a convolution 
of the observable's vector of binary strings (^f ) with the 
automaton's single cell images (j^^"^) and {\^^^)' 



\ tx^z ^z^z 
(tx^x *^x)(2^) 

itx^z*ex)ix) 




^ The definition of a scalar product is possible, because there are 
always only finitely many non-zero entries. 



The transformation ^ to the Laurent polynomials has 
the nice property of turning convolutions into multipli- 
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cations: 

y 

A;^a;+y l— — y 

= i-v- 

Because of this property, we refer to ([5]) as an algebraic 
Fourier transform. We get 

(ix\^)\ ^ ( ix^xin) iz^xiu) \ (ex{u)\ 
W\n)J \ ix-.ziu) iz^z{u) ) \eziu)) 

/ix^x (u) ■ i'x (u) + iz^x (u) ■ i'z (^)\ 
V^x^z{u)-ix (") + tz^z (m) • 4 {u))' 

In the following, we will omit the hat and the variable 
u for the sake of a short notation. Furthermore, we will 
replace £,x by £,z by and tx^x by t_n etc. to be 
consistent with the notation introduced in We then 
have 




Our example glider CQCA now looks like this: 

\1 u + u I 

We have already seen that the CQCAs have to ful- 
fill certain conditions, namely the local rule has to be a 
translation-invariant homomorphism which maps Pauli 
products to Pauli products. Furthermore, it has to obey 
commutation relations with its translates. We have to 
translate these conditions to the polynomial picture. The 
matrix we use does not have any dependence on the po- 
sition on the chain, so translation invariance is already 
included in this formulation. The commutation relations 
are encoded in the symplectic form, as we can see from 
(l3|). For T to conserve the commutation relations, the 
corresponding classical automaton has to conserve the 
symplectic form. This is why the classical CAs that corre- 
spond to CQCAs are called symplectic cellular automata 
(SCAs). In [l[ it was proven that a 2 x 2 matrix t with 
Laurent-polynomial entries is a SCA if and only if it ful- 
fills the following conditions: 

1. det(t) ^ aeZ; 

2. all entries are symmetric polynomials centered 
around the same (but arbitrary) lattice point a; 



3. the entries ti^, t2j of both column vectors, which 
are the pictures of (1,0) and (0, 1), are coprime. 

Furthermore, it was shown that to every CQCA there 
exists a SCA and an appropriate translation-invariant 
phase function A(^) such that 

T[w(C)] = A(Ow(te), (8) 

A(^ + 1]) = A(e)A(?/)(-l)«+''--(*«)+(*'')- 

and |A(^)| = 1 hold. Additionally, A(^) is uniquely 
determined for all ^ by the choice of A on one site. On 
the other hand, we can find CQCAs for any given SCA 
by adding a phase function. Thus, CQCAs and SCAs are 
equivalent up to a phase. We will therefore only refer to 
CQCAs, even if we talk about the corresponding SCAs. 

The last condition, the homomorphism property, is 
automatically fulfilled because the choice of the phase 
function and the conservation of the symplectic form 
a{t^,tri) = (t(^, 77). The multiplication of Pauli matri- 
ces is mapped to the addition (modulus two) of phase- 
space vectors. As our matrices are linear transforma- 
tions they obey t(^ 4- 77) = -|- trj. This translates to 
T([w(0w(,7)] - T[w(C)]T[w(7?)] via 

T[w(Ow(77)] 

= (-l)-''+«-T[w(e + r7)] 

= A(e + 77)w(t(e + r7))(-l)-''+«- 

= (-l)(*")+(*«)-w(tC)w(t77)A(e + ?7)(-l)-''+«- 

.(_!)?+')- -(t«)+(ti7)-(_l)-'?+«-+(tr))+(t5)_ 

= T[w(0]T[w(?7)](-l)'^(«'")~'^(*«'*'') 
= T[w(0]T[w(,7)]. 

A very simple CQCA is the shift on the lattice, which 
has the matrix It obviously commutes with all other 
CQCAs. Thus, we can multiply a CQCA with determi- 
nant u^"' which has entries that are centered around a by 
a shift by —a sites to obtain a centered CQCA with de- 
terminant 1 . From now on we will only consider centered 
CQCAs. 

A nice property of these centered CQCAs, which we 
will need to prove that the entanglement generation is 
linear, is that their matrices form a multiplicative group. 
Multiplying the matrices means concatenating the CQ- 
CAs. 



2. Classes of CQCAs 

CQCAs show a variety of time evolutions that can be 
roughly grouped into three classes. The first class shows 
periodic behavior, the second class consist of CQCAs that 
have glider observablcs, which just move on the lattice 
as shown in Figure [21 and the last case generates fractal 
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Figure 3: Time evolution of the fractal CQCA Tp with 
trtF ~ u"^ + 1 + u. Time increases upwards. The differ- 
ent colors mark the different Pauli matrices. 



space-time evolutions as shown in Figure [3l The class 
of a CQCA is determined by the trace of its matrix. If 
the trace is a constant trt = c, c S Z2, the automaton 
is periodic. If it is of the form trt = u^" + u" , n e N, 
the CQCA has gliders that move n steps on the lattice 
every time step. All other cases show fractal behavior 
0. We have already seen an example of a glider CQCA. 
The following CQCA exhibits fractal behavior: 



1 



1 



u 1 




(9) 



Its trace is trti? = u^^ + 1 + u. The time evolution that 
is shown in Figure [3] is very similar for all CQCAs with 
the same trace. This is observed also for other classes 
and subclasses of CQCAs. Mathematically, however, in 
general there is no notion of equivalence known for CQ- 
CAs with the same trace. An exception are CQCAs with 
trace u^^ + u, which can be shown to be all equivalent 
in the sense that they can all be transformed into the 
standard ghder CQCA ^ via 



(10) 



where b is a CQCA 



III. PURE TRANSLATION-INVARIANT 
STABILIZER STATES 

A. Definition 



Example III.l. The stabilizer group with the gener- 
ators § = {X ® X, Z ® Z} stabilizes the Bell state 
^ = l/-\/2(|l,l) + |0,0)). To check this, we just have 
to apply the stabilizer generators to the state (we omit 
the normalization): 

{X(g)X)i; = {X(E)X){\1)(E}\1)) + {X(E}X){\0)(E)\0)) 

= |0,0) + |1,1) 
iZ(E)Z)4> = |1)) + (Z®Z)(|O)0|O)) 

= (-|1))®(-|1)) + |0)®|0) 
= (-1)2|1,1) + |0,0) =v 



For infinitely many qubits we can't just write down the 
state, so we use an abstract definition. A translation- 
invariant stabilizer state is defined by a translation-in- 
variant set of operators § = {w(ra.^), x S Z}, where 
is the lattice translation by x sites. In [l[ it was proven 
that such a set defines a pure translation-invariant sta- 
bilizer state if and only if ^ is reflection invariant and 
the Laurent polynomials of ^ have no common divisors: 
gcd(^+,^_) = 1. For the polynomials, this means that 

1. ^ is of odd length, because reflection-invariant Lau- 
rent polynomials of even length are always divisible 
by 1 -I- M. We will write I ~ 2n+ 1. 



2. ^(0) 7^ (g); thus, the center element is not 1. 
^ would have the divisor + u. 



Else 



3. At least two different types of Pauli matrices (both 
different from the identity) have to occur (e.g. X 
and Y). Otherwise, ^+ = or ^_ = or ^_|- = 
each case implying common divisors. 

An example is Si ~ Xi-i ® Zi ® Xi+i. The condi- 
tion that the polynomials of ^ are coprime is also a 
condition for the column vectors of a CQCA matrix. 
This indicates a close connection between CQCAs and 
transla-tion-invariant stabilizer states. Indeed, CQCAs 
map pure translation-invariant stabilizer states onto each 
other. Furthermore, every translation-invariant stabilizer 
state can be generated by a single step of a CQCA from 
the "all-spins-up" state, which is a stabilizer state with 
stabilizer generators § ~ {w(tj;(0, 1)) = Z^, x G Z}. 
Thus, we can study the entanglement generation proper- 
ties of CQCAs acting on pure translation-invariant sta- 
bilizer states by just studying the entanglement in these 
states. 



A stabilizer state w is a common eigcnstatc to a group 
of commuting operators S = ({5^}). This means that 
u) o S = u}, yS & S. The group is generated by the 
set of generators § = {Si} by multiplication. It therefore 
suffices to check the stabilizer condition uj o S = uj for 
the generators 5*^. For finitely many qubits, this can be 
easily understood with the following example: 



B. Entanglement 

Entanglement in the stabilizer formalism was first 
studied in It was found that the entropy of en- 

tanglement with respect to a bipartite cut of some sta- 
bilizer state is exactly the number of generators of the 
correlation subgroup. The correlation subgroup Sab for 
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a bipartite stabilizer state is generated by all minimal 
stabilizer generators that have support on both parts 
of the system. The minimal stabilizer generators are 
those stabilizer operators which generate the stabilizer 
group (through multiplication) while having the small- 
est support of all such sets of operators. There are in 
general a lot of such sets, so the choice is not unique. 
For example, the stabilizer generators Z (E) 1 and t (E) Z 
would stabilize the product state "0 = It) ® It)- Obvi- 
ously the combination Z (E) Z and t E) Z would do the 
same, but with one generator with larger support. In 
this case the first set would be minimal, while the second 
wouldn't. In the case of translation-invariant pure stabi- 
lizer states we don't have this issue. The set of generators 
is always translation invariant and fulfills the conditions 
introduced in Section IIII Al Assume that a given gen- 
erator of a translation-invariant pure stabilizer state is 
not minimal. Then it is composed of at least two stabi- 
lizer generators which also have to fulfill the conditions 
from Section IIII Al In particular, all of the generators 
have to be the same. This implies that the polynomi- 
als of the original non-minimal generator have common 
divisors. But this is not possible, because it is required 
that the polynomials are coprime. Thus the generators 
of translation-invariant pure stabilizer states are always 
minimal. 

Unfortunately the proof in [l^ relies heavily on the 
fact that only finite systems are considered. We use a 
different approach and come to essentially the same result 
that holds for translation-invariant stabilizer states on 
infinite chains. 



1. The Bipartite Case 

First we will investigate the case of a bipartite splitting 
of the chain. We have two parties, say Alice and Bob, 
where Alice controls the part and Bob the part 3§ of 
the chain. This is shown in Figure |4l Let us first define 



iHuuy 



□□□□□ 



Figure 4: The spin chain is cut into two halfchains jz/ and , 
We study the entanglement between these halfchains. 



In this case we have the following theorem: 

Theorem III. 3. A pure translation-invariant stabilizer 
state of stabilizer generator length 2n+l entangles n qubit 
pairs maximally with respect to any bipartite cut. 

Proof. Here we will only present the idea of the proof. 
The technical parts can be found in the appendix of 

Obviously, all bipartite cuts are equivalent as the state 
is translation- invariant. Thus, we can look at any partic- 
ular cut to prove the general result. We have a stabilizer 
generator centered around each site. Unless our stabilizer 
generators are single site operators (n = 0), the cut will 
always leave several stabilizer generators cut into parts 
on both systems ^ and As one can see in Figure [5] 
there will be 2n stabilizer generators affected. These op- 
erators generate the correlation subgroup Sab- Thus, 
the correlation subgroup has 2n generators and there- 
fore dimension |5yiB| = 2n. The interesting fact about 
these stabilizer generators is that, despite commuting as 
a whole, their restrictions to ox don't necessarily 
commute. Now we try to find commuting pairs of anti- 
commuting Pauli products in the restriction of Sab to 
or ^ by multiplying stabilizer generators from the corre- 
lation subgroup (We check for the anti-commuting parts 
only on one side, but carry out the multiplication on both 
sides to preserve the stabilized state). We know from the 
theory of quantum error correction codes that each such 
pair on .s^ encodes one qubit [§]. Because we carried 
out the multiplication on both sides, the corresponding 
parts of the operators on SB fulfill the same commutation 
relations and therefore also encode qubits. The pairs of 
operators on each halfchain behave like pairs of X and Z\ 
thus, we call them X and Z. We only required them to 
anti-commute, so we did not fix which of them is the X 
and which the Z. Thus, we can choose this and we choose 
it in such a way that the corresponding operators on each 
side are either bot Z or both X. Thus, we have a pair of 
stabilizer generators that reads E) Xag and Z,^ E) Zag. 
But, as seen in Example IIII.H ^ = {X (E X,Z (E Z'] en- 
codes a Bell state, which is maximally entangled. Thus, 
each pair of anticommuting pairs stabilizes a maximally 
entangled (logical) qubit pair. 

What is left to show is that we can always find n such 
pairs. The proof is rather lengthy and technical. It is 
carried out in Q and based on methods from quantum 
error correction codes to directly construct the pairs. □ 



bipartite entanglement for stabilizer states. 

Definition III. 2. The entanglement E{uj^) of a 
translation-invariant stabilizer state oj^ in a bipartite set- 
ting is the number of maximally entangled qubit pairs with 
respect to any bipartite cut. These qubit pairs are logical 
ones, that is they are each localized on several physical 
qubits. By local operations on each part of the chain one 
could localize the logical qubits onto one physical qubit 
each. 



C. The Tripartite Case 

In this setting we cut the chain into three parts, one 
middle part of length L and two infinite ends, as shown 
in Figure [6l We now want to calculate the entanglement 
between the finite part and the two infinite parts. To do 
this calculation, we use the same method as above, and 
arrive at the following theorem. 
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Figure 5: Cut stabilizer generators: it is apparent that for 
operators of length i = 2n + 1, 2n of them have support on 
both and SS. 



□ □□□ 



Figure 6: A finite region .g/ of 4 spins is cut out of the chain 
leaving two infinite ends and . 



Theorem III. 4. Given a pure translation-invariant sta- 
bilizer state uj^ of stabilizer generator length 2n + 1, a 
region of length L shares 2n maximally entangled qubit 
pairs with the rest of the chain if 2n < L and L qubits 
pairs if 2n > L. 

Proof. The proof works exactly as in the bipartite case. 
In the case 2n < L the cut stabilizers are only cut on one 
side. But all stabilizers that are cut on the left-hand side 
commute with those cut on the right-hand side. Thus, 
we have two independent cuts of the bipartite case and 
therefore 2n pairs of maximally entangled qubits. In the 
case 2n > L some stabilizers arc cut on both sides. We 
use the same technique to produce the mutually com- 
muting anti-commuting pairs which encode the qubits as 
in Theorem IIII.3I We always find L pairs of maximally 
entangled qubits. □ 



IV. ENTANGLEMENT GENERATION 

Now we come to the generation of entanglement 
through CQCA action. As we have seen in the previous 
section, the bipartite entanglement of a pure translation- 
invariant stabilizer state depends linearly on the length 
of the generators of the stabilizer group. So, it suggests 
itself to study the evolution of the length of the stabilizer 
generators under CQCA action. First, let us define the 
asymptotic entanglement generation rate from stabilizer 
states. 

Definition IV. 1. The asymptotic entanglement genera- 



tion rate from stabilizer states for CQCAs is defined as 

(11) 



AE . . , 1 , 



where E{LU^,t) is the bipartite entanglement at time t. 

Lemma IV. 2. The length 2n + 1 of the minimal stabi- 
lizer generators grows asymptotically with 

= lim 7(2n(^5,t) + 1) = dg(trt) (12) 
Ar t— i-oo t 

for any centered CQCA T and any translation-invariant 
pure stabilizer state lo^ . 



Proof. Wc know that CQCAs map pure translation-in- 
variant stabilizer states onto pure translation-invariant 
stabilizer states. The image of a state with stabilizer 
generators § = {'w{tx^), x £ Z} under the action of 
t steps of a CQCA T is a state with stabilizer gen- 
erators §* = {w(r3;t*^), X e Z}. Furthermore, we 
know that any stabilizer state can be generated from 
the "all-spins-up" state by a CQCA b. So we have 
S* = {w(r:j;t*b(0, 1)), X e Z}. The length of the sta- 
bilizer generator is determined by the highest order of 
the stabilizer generator polynomials, dg(^). Namely, the 
stabilizer generator is of length 2 • dg(.^) -|- 1. So, we have 
to calculate dg(t*^) = dg(t*b(°)). 

For an arbitrary product of CQCAs HiLi '^^ '^^^ 
define the series (a/)i</<fc = dg(J|'^^ti). It is subaddi- 
tive, i.e. a„_|_m < a„ -|- a^, because the concatenation of 
CQCAs is essentially the multiplication and addition of 
polynomials, which is subadditive in the exponents. For 
subadditive series a„ Fekete's lemma [lH states that the 
limit lim„_>.oo ^ exists. In our case the series is always 
positive, so the limit is positive and finite. An easy way 
to determine the limit is to take a subseries, which of 
course has the same limit. The subseries of the t = 2'^th 
{k £ N) steps is a good candidate, because we can make 
use of the Cayley-Hamilton theorem to obtain 



t^' =t(trt)2'=-i-f l^(trt)2'=-2\ 
Furthermore, we have 

dg(t2^bO) ^ dg(tb(°)(trt)^^-i+b(°)^(trtf^-2') 
= a ■ c{k) + b ■ d{k), 

with a = tbC^), b = b(°), c(fc) = (trt)2'-i, d{k) = 

Eli(trt)2^~-2\ 

Let us first assume that for some ko we have 
dg(t^ b(i)) ^ •ig(6). We start by determining a recur- 
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sion relation for c(fc) and d{k). 



c(fc + 1) 
d{k + 1) 



(trt)"'"-^ = (trt)' 
c(fc)(trt)2' 



^(trt)^ 



fc+i 



1=1 



(trtf 



= d{k){tTtf' + 1. 

Now we are able to calculate the limit: 
= lim ^dg{a ■ c{k) + b ■ d{k)) 



lim 



2k- 



-dg(a.c(fc-l)(trt) 

fc-i-oo Z 

+&• l)(trt)2'"' +b) 

lim ldg((trt)2'"r(fc-l)) 

fc— >CX3 Z 

'2 



i™, ^ ( dg(trt)- + dg(r(fc - 1)) 



1™. ( E ^dg(trt) 



^dg(r(fco)) 



= dg(trt). 

In the third step we used that dg(6) < dg(t2'' b(°)), Vfc > 

In the second case we have dg(t2'"b(°)) < dg(6), Vfc. 
Thus, dg(t2'b(°)) is bounded by dg(fe). Therefore, n{t, £,) 
is bounded and = 0. But n{t, ^) bounded also implies 
t periodic and therefore dg(trt) = 0. Thus, we have 
^ = = dg(trt). This completes the proof. □ 



Now we only have to put Theorem IIIL3I together with 
Lemma riV.2l to prove the following theorem: 

Theorem IV. 3. The asymptotic bipartite entanglement 
generation rate (maximally entangled qubit pairs per 
step) of a general centered CQCA T is the degree of its 
trace polynomial, dg(trt). 



Proof. By Theorem lIII.3[ every translation-invariant sta- 
bilizer state with stabilizer generator length 2n + I en- 
tangles n qubit pairs maximally with respect to any bi- 
partite cut. Lemma llV.21 shows that under the action of 
a CQCA T the length of the stabilizer generators grows 
asymptotically with dg(trt) qubit pairs per step. This 
implies the proposition. □ 

This means that starting from the "all-spins-up" stabi- 
lizer product state, the entanglement grows linearly with 



dg(trt) under the action of a CQCA T, because b = 1. 
Starting from an arbitrary translation-invariant pure sta- 
bilizer state, the entanglement might decrease in the first 
k steps, e.g. if t'^ = b~^, but then starts to increase lin- 
early with dg(trt). This behavior is shown in Figure [T] 



40 



30 



20 



Tg, tr — Ti ^ + w 
2-step glider CQCA, tr = u^' + 

Tf, tr = + 1 + u 
periodic CQCA, tr = 



10 

t 



Figure 7: Entanglement generation for the stabilizer state 
with w(5) = y®X®Fina bipartite setting with differ- 
ent CQCAs. One can see that entanglement can also be de- 
stroyed, but grows asymptotically linearly with the number 
of time steps t. The rate is given by the degree of the trace 
of the CQCA matrix. 



In the tripartite setting, again starting from the "all- 
spins-up" state, the entanglement between the finite re- 
gion and the rest of the chain grows under the action of 
a CQCA T with 2 • dg(trt) until it reaches L. Then it 
remains constant. If we start with a general translation- 
invariant stabilizer state, again the entanglement might 
decrease at first. After some time it starts increasing and 
reaches L, where it remains if the CQCA is not periodic. 
The entanglement generation is twice as fast as for the 
bipartite setting, because we have two cuts. Results are 
shown in Figure [8l 

In both cases, the asymptotic entanglement generation 
only depends on the degree of the trace of the automa- 
ton. A CQCA with gliders generates entanglement just 
as fast as a fractal automaton if the degree of the trace 
is the same. Of course, periodic automata destroy all en- 
tanglement they just generated in the next steps. This 
corresponds to the fact that their trace is a constant. The 
entanglement generation rate is however not directly gov- 
erned by the neighborhood of the CQCA. Of course, the 
size of the neighborhood bounds the possible rate of en- 
tanglement generation from above, but even automata 
with a huge neighborhood can be periodic and thus gen- 
erate no entanglement at all. 

It is worth mentioning that CQCAs saturate the bound 
on the entanglement generation rate for translation- 
invariant operations acting on translation-invariant 
states derived in This means that there is no 
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Tg, tr — u ^ + u T 

: 2-step glider CQCA, tr ^ + a 

Tf, tr — + 1 + u x 

periodic CQCA, tr — q 

30 - A ■5' ^- ^ ^- — -*"-"«—-«> M -W K - 
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^ X 

□ ^ a _<si" X E3 El H □ ■ 
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5 10 15 20 

t 

Figure 8: Evolution of entanglement of a subchain of 30 
consecutive spins for an initial stabilizer state with w(^) = 
Y®X®X®X®X®X®Y for different CQCAs. The 
entanglement first grows as in the bipartite case, but then 
saturates at 30 qubit pairs. 



ends of the spin chain. However, for special cases the 
cut images still fulfill the commutation relations. Fortu- 
nately, the much used glider CQCA is of this type. So, 
in this case the only influence from the ends of the chain 
is that the outermost sites miss the influence from one 
neighbor. As we still have an automorphism and thus a 
reversible operation on the whole chain, no information 
can be lost at the ends of the chain. So, the ends have 
to be reflective. We can observe that in the case of an 
incoming glider in the left part of Figure [51 This leads to 
the fact that observables are mirrored by the chain. For 
single-site Pauli matrices, this can be seen in the right- 
hand part of Figure [9l Arbitrary observables are sums of 
tensor products of single-site Pauli matrices. Because the 
single-site Pauli matrices are mirrored, the tensor prod- 
ucts will be mirrored, too. Thus, also sums arc mirrored 
and therefore any observable is mirrored on the chain. 
This property and the spreading of observables in this 
process is used in applications of CQCAs like 0, ll]. 

VI. APPLICATIONS 



translation-invariant operation that generates more en- 
tanglement per step, while having the same size of neigh- 
borhood, than a CQCA whose trace has maximal degree 
with respect to the neighborhood. 



V. NOTES ON FINITE CHAINS 

In our whole analysis we used an infinite spin chain 
to obtain translation invariance, which would be broken 
by the ends of the spin chain. For finite chains we gen- 
erally have to take the effects of the ends into account. 
But when we deal with observables which are localized 
far away from the ends of the chain and time periods T 
which are much smaller than the length of the chain, for 
locality reasons the time evolution of the observables has 
to be the same as for the infinite chain. Information from 
the ends of the chain can only travel with the same finite 
speed, as all other information, so e.g. a one-site observ- 
able localized in the middle of a chain of length L = 21 + 1 
can only interfere with influences from the ends after 1/2 
time steps. 

Admittedly the applications of CQCAs make heavy use 
of the effects occurring at the ends of the spin chains. 
Thus, our above statement is not applicable. To take ef- 
fects at the ends of the spin chains into account, we have 
to deal with the fact that they break the translation in- 
variance. The one-site images of the sites at the ends of 
the chain whose neighborhoods would reach over one end 
of the chain, if they were the same as on the rest of the 
chain, have to be adapted. Like in the case of stabilizer 
states above, in general the cut images don't fulfill the 
necessary commutation relations any more. There is no 
general theory yet of how to adapt the CQCAs at the 



CQCAs are used in different quantum computational 
schemes. Of course they are not capable of doing any 
computations efficiently; a classical computer couldn't 
do efficiently as well, because they are essentially classical 
cellular automata. But, nevertheless they can be of great 
use when accompanied by non-Clifford operations. This 
is used in different approaches. 

A. Measurement Based Quantum Computation 

The most famous model of quantum computation in- 
volving a CQCA (indirectly) is the idea of measurement 
based quantum computation or the "one way quantum 
computer" by Raussendorf and Briegel @ . The resource 
state used to perform the quantum computation by suc- 
cessively measuring the single qubits is a so-called cluster 
state on a 2-dimensional (finite) lattice of qubits. It can 
be generated by a CQCA, namely a 2-dimensional version 
of the ghder CQCA Thus, one time-step of a CQCA 
together with measurements of single qubits suffices for 
universal quantum computation. 

B. Raussendorf 's Scheme of Translation Invariant 
Quantum Computation 

In this scheme of translation-invariant quantum com- 
putation 0, again the glider CQCA is used, in this case 
the one-dimensional version. The property of generat- 
ing patterns from single-site observables, thus spreading 
them over the spin-chain (a finite chain is used here), 
is used to immunize observables against special global 
transformations. The time steps in which certain ob- 
servables arc immune against these operations depend 
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Figure 9: Time evolution of the glider CQCA on a finite chain of 7 spins. In the left-hand illustration we can see that the 
gliders are reflected at the ends of the chain. In the other two illustrations we see that on a finite chain the glider CQCA 
mirrors the position of single-site Pauli matrices. As all Pauli products can be decomposed into single-site Pauli matrices, this 
holds for all observables. (All observables can be decomposed into a sum of Pauli products, which are each mirrored in the 
same number of time steps, so the whole observable is mirrored.) 



on their initial position on the chain. Thus, temporal 
control (when to apply the global gates) can be turned 
into spatial control and any quantum operation can be 
conducted on the system. An example of such behavior 
would be the effect of a global Pauli Y operation (a local 
Y on each site). It gives a sign on each X and Z tensor 
factor. So, in our example shown in Figure [9] the Z-2 
would gain a phase of —1 when the Y is applied in one 
of the steps {0,1,6,7,8}, while it would gain no phase 
in the other steps. The X_i gains the phase in steps 
{0, 1, 2, 3, 6, 7, 8}. The steps where a phase is gained be- 
long to the contraction resp. expansion of the observables, 
while the steps where no phase is gained belong to the 
transmission of the expanded observable over the chain. 
So, if we e.g. run the automaton for 2L + 2 = 16 time- 
steps and apply a global Y in step 3, X-i gains a phase 
but Z-2 does not. Applying Y in steps {1,3,6,7,} Z-2 
gets the phase. So, we turned temporal control into spa- 
tial control. To achieve universal quantum computation, 
we use arbitrary translation-invariant local rotations in- 
stead of the global Y. 

C. The Quantum Computational Scheme by 
Fitzsimons and Twamley 

In Q the same property of the glider CQCA is used. 
However in this case the non-Clifford operations are not 
translation-invariant, but only conducted at the ends of 
the chain. Separate control of the ends of the chain is 
justified by the fact that due to the missing neighbor the 
physical properties of the systems at the end differ from 
those in the middle of the chain. The CQCA is used to 
transport qubits to the end of the chain, which can then 
be manipulated. Two-qubit gates are achieved by first 



decoupling one spin (at the end of the chain) from the 
chain and then transporting the other one to its neigh- 
boring position. Then the gate is applied and the qubits 
are transported back. This scheme was experimentally 
realized in a NMR-system [T2j . 



Conclusion and Outlook 

We have introduced CQCAs and their classical descrip- 
tion. It was shown that CQCAs acting on translation- 
invariant stabilizer states generate entanglement with the 
highest possible rate for translation-invariant operations. 
The rate only depends on the trace of the CQCA's ma- 
trix. It is independent of the class of the CQCA. We 
furthermore commented on some applications of CQCAs 
in quantum computational schemes. These schemes use 
finite chains of qubits, while our analysis is based on in- 
finite spin-chains. Thus, a future task is, to complete the 
theory of CQCAs to include finite spin chains. This was 
already done for periodic boundary conditions [ij, but 
the case of non-periodic boundaries is still an open ques- 
tion. Additionally, the investigation of invariant states 
and convergence of states under CQCA action which 
was started in Q will be continued for finite chains to 
strengthen the connection to the applications of CQCAs. 
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